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<N ; Abstract 

C^ , Applying effective Lagrangian method and on-shell scheme, we analyze the electroweak correc- 

tions to the rare decay 6 — )■ s + 7 from some special two loop diagrams in which a closed heavy 
CN . fermion loop is attached to the virtual charged gauge bosons or Higgs. At the decoupling limit 

(— I \ where the virtual fermions in inner loop are much heavier than the electroweak scale, we verify the 

final results satisfying the decoupling theorem explicitly when the interactions among Higgs and 






heavy fermions do not contain the nondecoupling couplings. Adopting the universal assumptions 
on the relevant couplings and mass spectrum of new physics, we find that the relative corrections 



^ ' from those two loop diagrams to the SM theoretical prediction on the branching ratio of S — t- X 7 

"^ ' can reach 5% as the energy scale of new physics Aj^p = 200 GeV. 
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I. INTRODUCTION 

The rare B decays serve as a good test for new physics beyond the standard model 
(SM) since they are not seriously affected by the uncertainties originating from long distance 
effects. The forthcoming and running B factories will make more precise measurements on 
the rare i?-decay processes, and those measurements should set more strict constraints on 
the new physics beyond SM. The main purpose of investigating 5-decay, especially the rare 
decay modes, is to search for traces of new physics and determines its parameter space. 

The measurements of the branching ratios at CLEO, ALEPH and BELLE |JL] give the 
combined result 

BR{B -^ X 7) = (3.52 ± 0.23 ± 0.09) x 10"^ , (1) 

which agrees with the next-to-next-to-leading order (NNLO) standard model (SM) predic- 
tion y 

BR{B -> X^7) = (3.15 ± 0.23) x 10"^ . (2) 

Good agreement between the experiment and the theoretical prediction of the SM implies 
that the new physics scale should lie well above the electroweak (EW) scale A^,^. The 
systematic analysis of new physics corrections to 5 — )■ X^7 up to two- loop order can help 
us understanding where the new physics scale sets in, and the distribution of new physical 
particle masses around this scale. In principle, the two-loop corrections can be large when 
some additional parameters are involved at this perturbation order besides the parameters 
appearing in one loop results. In other words, including the two-loop contributions one 
can obtain a more exact constraint on the new physics parameter space from the present 
experimental results. 

Though the QCD corrections to the rare B decays are discussed extensively in literature, 
the pure two-loop EW corrections to the branching ratio of 6 —t- 57 are less advanced 
because of the well known difficulty in calculation. Strumia has evaluated the two-loop 
EW corrections to 6 — )■ S7 from the top quark using heavy mass expansion in gaugeless 
limit of the SM ^ . At the limit of large tan (3 in supersymmetry, Ref . [^ analyzes the two 
loop corrections to the branching ratio of S — > X^7 from the virtual charged Higgs and 
gluino-squark sector. 

Employing the effective Lagrangian method and on-shell scheme, we present the correc- 
tions to the branch ratio of -B — ;■ X^7 from some special diagrams in which a closed heavy 



ferinion loop is attached to the virtual charged gauge bosons or Higgs here. The effective 
Lagrangian method can yield one loop EW corrections to the effective Lagrangian of 6 — )■ S7 
exactly in the SM and beyond, and has been adopted to calculate the two loop supersym- 
metric corrections for the branching ratio of 6 — ;■ 57 |5|, neutron EDM [6] and lepton MDMs 



and EDMs 



Qls 



8|. In concrete calculation, we assume that all external quarks and photon 
are off-shell, then expand the amplitude of corresponding triangle diagrams according to 
the external momenta of quarks and photon. Using loop momentum translating invariance, 
we formulate the sum of amplitude from those triangle diagrams corresponding to same 
self energy in the form which explicitly satisfies the Ward identity required by the QED 
gauge symmetry, then get all dimension 6 operators together with their coefficients. After 
the equations of motion are applied to external quarks, higher dimensional operators, such 
as dimension 8 operators, also contribute to the branching ratio of i? — )■ X^7 in principle. 
However, the contributions of dimension 8 operators contain the additional suppression fac- 
tor IV? / h^ comparing with that of dimension 6 operators, where m^^ is the mass of bottom 
quark. Setting A^,^ ~ lOOGeV, one obtains easily that this suppression factor is about 10^^ 
for the h — )■ 57. Under current experimental precision, it implies that the contributions of 
all higher dimension operators {D > 8) can be neglected safely. 

We adopt the naive dimensional regularization with the anticommuting 7g scheme, where 

there is no distinction between the first 4 dimensions and the remaining D — A dimensions. 

Since the bare effective Lagrangian contains the ultraviolet divergence which is induced 

w divergent subdiagrams, we give the renormalized results in the on-mass-shell scheme 



9|. Additional, we adopt the nonlinear R^ gauge with .^ = 1 for simplification 10|]. This 
special gauge-fixing term guarantees explicit electromagnetic gauge invariance throughout 
the calculation, not just at the end because the choice of gauge- fixing term eliminates the 
jW^G^ vertex in the Lagrangian. 

This paper is composed of the sections as follows. In section [TTl we introduce the effective 
Lagrangian method and our notations. We will demonstrate how to obtain the identities 
among two loop integrals from the loop momentum translating invariance through an ex- 
ample, then obtain the corrections from the relevant diagrams to the effective Lagrangian 
of 6 — 7- 57. Section [Till is devoted to the numerical discussion under universal assumptions 
on the parameters of new physics. In section \1V\ we give our conclusion. Some tedious 
formulae are collected in the appendices. 



II. THE WILSON COEFFICIENTS FROM THE TWO-LOOP DIAGRAMS 



In this section, we derive the relevant Wilson coefficients for the partonic decay 6 — )■ S7 

including two-loop EW corrections. In a conventional form, the effective Hamilton is written 

as 
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where V is the CKM matrix and G^, = 1.16639 x 10 ^ GeV ^ is the 4-fermion coupling. 



The definitions of those dimension six operators are 11 1 
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where F^^ and G ^ = G"'^jr°' are the field strengths of the photon and gluon respectively, 
and T" (a = 1, ■ ■ ■ , 8) are SU{?>)^ generators. In addition, e and g^ represent the EW and 
strong couplings respectively. 
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FIG. 1: The relating two-loop diagrams in which a closed heavy fermion loop is attached to virtual 
W^ bosons or G^ {H^), where a real photon or gluon is attached in all possible way. 

After expanding the amplitude of corresponding triangle diagrams, we extract the Wilson 
coefficients of operators in Eq. (jl]) which are formulated in the linear combinations of one and 
two loop vacuum integrals in momentum space, then obtain the corrections to the branching 
ratio of B -^ ^sl- Taking those diagrams in which a closed heavy fermion loop is inserted 
into the propagator of charged gauge boson as an example, we show in detail how to obtain 
the Wilson coefficients in effective Lagrangian. 

A. The corrections from the diagrams where a closed heavy fermion loop is in- 
serted into the self energy of W^ gauge boson 

In order to get the amplitude of the diagrams in Fig{T]^a), one can write the renormal- 
izable interaction among the charged EW gauge boson W^ and the heavy fermions F^^^ in 
a more universal form as 



c 



-1^-''^F„7m(C>- + Cu^+)Fp + h.c. 



(5) 



where the concrete expressions of (^'^ depend on the models employed in our calculation. 
The conservation of electric charge requires Q^ — Qa = 1, where Qa,/3 denote the electric 
charge of the heavy fermions Fq,^^ respectively. 

Applying Eq.(l5]), we write firstly the amplitude of those two loop diagrams in FigIT](a). 
For example, the amplitude for the diagram in which a real photon is attached to the virtual 
W^ boson (Figj2]) can be formulated as 
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FIG. 2: The triangle diagram in which the real photon is attached to W^ gauge boson. The 
amplitude is written in Eq([6]). 
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Here Aj^^, denotes the renormalization scale that can take any value in the range from the 
EW scale A^^ to the new physics scale Aj^p naturally, and we adopt the abbreviations: 



cos 9„ 



sin6'„ with 6^ denoting the Weinberg angle. Additionally, p, k are the 



incoming momenta of quark and photon fields, p is the Lorentz index of photon, separatively. 
Certainly, the amplitude does not depend on how to mark the momenta of virtual fields 
because of the translating invariance of loop momenta. 

It can be checked easily that the sum of amplitude for diagrams in Fig{T]^a) satisfies the 
Ward identity required by the QED gauge invariance 



kPA^^^Hp, k) = efsffi'^))(P + k)- sffi'^)) 
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where ^Wa)) denotes the sum of amplitudes for the triangle diagrams corresponding to the 
self energy in Fig{T]^a), as well as Sj^'^-''' denotes the amplitude of corresponding self energy 
diagram, respectively. 

According the external momenta of quarks and photon, we expand the amplitude in 
Eq.® as 
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since we only consider the corrections from dimension 6 operators, here P„„ = {qi—Tn^^){qi — 
miy{{q2 - qif - mlj{ql - m^ . 

Because the denominators of all terms in Eq.([H]) are invariant under the reversal gi -^ 
—qi,q2 — ^ —q2, those terms in odd powers of loop momenta can be abandoned, and the 
terms in even powers of loop momenta can be simplified by 
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and those similar formulae presented in Eq.(5) of Ref[5[, where the tensors are defined as 
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Summing over those indices which appear both as superscripts and subscripts simultane- 
ously, we derive all possible dimension 6 operators in the momentum space together with 
their coefficients which are expressed in the linear combinations of one and two loop vacuum 
integrals. In a similar way, one obtains the amplitude of other diagrams. Before integrating 
with the loop momenta, we apply the loop momentum translating invariance to formulate 
the sum of those amplitude in explicitly QED gauge invariant form, then extract the Wilson 
coefficients of those dimension 6 operators listed in Eq.(jl]). Actually, we can easily verify 
the equation 

d^qi d^q2 qi^, ^ , , 

{2n)D {2n)^ {q! - mlM - niMq, - q,y - ml) - " ^ ^ 

Performing an infinitesimal translation gi — > gi, q2 ^ 12 ^ o- with Op — t- (p = 0, 1, ■ ■ ■ , D), 

one can write the left-handed side of above equation as 
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In a similar way, other identities presented in Ref.|5[ can be derived. Using the expression 
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which is equivalent to the identity Eq.([T3l). Here, e = 2—D/2 with D denoting the dimension 
of space-time, A is a energy scale to define Xj = mf/h? and x^ = A^ /A^. Additionally, 
A^ = x"^ +y'^ + z'^ — 2xy — 2xz — 2'yz, and the concrete expression of \E'(x,?/, 2;) can be found in 
the appendix. Actually, the equation Eq. flT6|l provides a crosscheck of Eq. flT^ and Eq. flT5|) 
rather than a verification of Eq.( fT3l) . After applying those identities derived from loop 
momentum translating invariance, we formulate the sum of amplitude from those triangle 
diagrams corresponding to the self energy FigHJ^a) satisfying QED gauge invariance and 
CPT symmetry explicitly, and extract the Wilson coefficients of those operators in Eq.(|l]). 
Integrating over loop momenta, one gets the following terms in the effective Lagrangian: 
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where a^ = e^ /An and Q^ = —1/3, Q^ = 2/3 represent the charge of down- and up- 
type quarks, respectively. T^ = 1 when the heavy virtual fermions take part in the strong 
interaction, otherwise T'^ = 0. The functions ipi, 0, are defined as 
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Note that the result in Eq{T7] does not depend on the concrete choice of energy scale A, 
and the concrete expressions of Fi{x,y,z,u), g. {x,y) {i, j = 1, 2 ■ • •) can be found in 
appendix. 

The charged gauge boson self energy composed of a closed heavy fermion loop induces the 
ultraviolet divergence in the Wilson coefficients of effective Lagrangian, the unrenormalized 
W^ self energy is generally expressed as 

,2 

where the form factors A^^ 2 ^-^^ 3^2 only depend on the virtual field masses and renor- 
malization scale. Here, we omit those terms which are strongly suppressed at the limit of 
heavy virtual fermion masses. The corresponding counter terms are given as 
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The renormalized self energy is given by 
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For on-shell external gauge boson W , we have 
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where e(p) is the polarization vector of W^ gauge boson. Inserting Eq. ( TT9l) and Eq. (jJO 
into Eq. ( l22l) , we derive the counter terms for the W^ self energy in on-shell scheme as 
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To cancel the ultraviolet divergence and those dangerous terms violating the decouphng 
theorem completely, we should derive the counter term for the vertex '-fW~^W~ here since 
the corresponding coupling is not zero at tree level. In the nonlinear R^ gauge with ^ = 1, 
the counter term for the vertex 'jW~^W~ is 
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where ki {i = 1, 2, 3) denote the incoming momenta of W^ and photon, and /x, u, p denote 
the corresponding Lorentz indices respectively. 

We can verify that the sum of amplitude from counter diagrams satisfies the Ward identity 
required by the QED gauge invariance obviously. Accordingly, the effective Lagrangian from 
the counter term diagrams is written as 
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•^ T^^ ~r a;. 



1 



144 ^ ^7(3;^, a; J + ^02(a;^,,a;^^)^5(a;„,a;J + ^8(a;^,a;J 



Co 



, /^L*^/? I ^H*Ai Va; a; U/2 



Y^V'2(a;„,a;J 



+Y^^i,i (a^F,, 2:^^)^2(3;^, a; J + 03(a;^„ , x^^)V^i(x„, xj + ^^{x^.x^) 



( 1 1 



a'Fr, ' '^Ffl/r'6(,a;„; Xj j 



+04(3;^ , a;„ )V'5(a;„, xj + ^io(x„, xj 



O. 



1 X ~\~ X 

+fC^*C^ +C^*C^) — fe + (a;p +x,Uq](x ,x,) + — ^Mx ,x, 

V a/3 ^a/3 ^a/3 ^a/3y 24^^ ^ -Fa F^ / Y'O J V w 5 tj ' \^/^ ^ "' * 
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Q-zA^F^^^pMei^^^^t) + oM^F^^^EMd^^^^t) +Mx^,Xt) 



24 



O. 



'y^ap'-'aP ' ''a/3 '''alsJ^^Fa-'^Fij) 



1/2 



Y^M^^^^t) + -^QlA^F^^^FrJMx^^Xt) 



fa ' i-p ' 



+04(x^^ , Xp^ )V^5 (x„ , X J + V^io(a;„ , x, ) 



O. 



(25) 



Adding the counter terms to bare Lagrangian Eq. fITTj) , we cancel the ultraviolet divergence 
there. Under our approximation, the resulted effective Lagrangian is written as 



- rr \/2G„a X 

pejf _ V F e ^ Y*Y J 

V^a/3 ^a/3 ^afi ^a, 

V^c^ ^a/3 ^Q/3 ^a, 

xSafl^afl ^a/3 ^a, 

\(aL*/-L _i_/-R*/-E 

V "^ a/3 ^a/3 ^a, 

V^c^ ^a/3 ^Q/3 ^c 

V a/3 a/3 ^a/3 ^( 

V a^ a/3 ^a/3 ^< 



^L* /-L 



-R*/-R 



^a,S ^Q/3 ^Qf/3 ^Q/3 



i^i + Q.F, 



l-^w 5 -^ t 5 -^ Fa 1 "^Fg ) ^2 



^ 'iK-^vjl -^tl -^Fa ' •^Ffl/'-^2 



yXp^ Xp^ ) 



1/2 



F, + Q^F, 



l-^w 5 -^t 5 '^F„ ' -^Fn J ^2 



l-^Fa-^F^J -Tgl^X^, Xj, X^^ , X^^ jCg 



^2 + t:f. 



\Xp^ Xp^ j 



l-^w' -^t' "^Fa' -^F^ J^6 



F, + T^F, 



l-^w 5 -^t 5 -^Fq ' "^F^ J '-^6 



-t _ -Tq I X„, , X. , Xc' •Xc'^JL/, 



9 V^wJ -^t; -^Fc ' Fp , 



X 



^a'^Fp) -'a-^iol-^w! ^t1 ^Fa ' '^F^J^12 f + ' 



(26) 



which only depends on the masses of virtual fields. It should be clarified that the corrections 

to the coefficients of Oy ^2 '^o not depend on the concrete renormalization scheme adopted 

here since the relevant terms from bare Lagrangian do not contain the ultraviolet divergence. 

In the limit z <^ x,y, the function $(x, y, z) can be approximated in powers of z as 

z z 

$(x, y, z) = ^o{x, y) + zipi{x, y) + ^<^2(a;, y) + ^</'3(a;, y) 



1 Q 

+2z(\nz - l)'K,{x,y) + 2z'(^ - ^Tx,{x,y) 
., /In^ 11 \ , 



(27) 



with 



7r,(x,?/) 
71-2 (a;,!/) 

T^i{x,y) 



X + y 2xy 



In 



{x — yY {x — yY x' 
1 12x?/ 



6xy(x + 1/) jj^ ^ 



{x-yY {x-yY (x-y) 



X 



(28) 
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and the concrete expressions of function ipi{x,y) {i = 0, 1,2,3) can be found in appendix. 
Using the asymptotic expressions in Eq.( l27|) . we derive the leading contributions contained 
in Eqj26] under the assumption m^ = rrip = uip ^ m„ : 






^slQ, 



X 



1 - 3Qff d^g^, _ 2 - 3Qfs dg,, _ 1 d^g^, _ 1 ^^^3,^ 

8 (9x2 8 Q^ X44 Q^4 43 Q^3 

w w w w 

1 d^g^, 1 d^g,, 29d'^g,,, 11 dg,. 



144 dx^ 



12 dx^ 72 (9x2 






12 (9x 



J l^;^, x^j 



O. 



[a^w.a^JCa 






1 (9^^4i 1 d^g,, 152^ 1 5ft n. ^ 



144 (9x4 



+ 



144 (9x4 ig ^3,3 

w w 

1 d^0s,i 5 g^ft^i 
24 9x2 



4 9x2 ]^g ^2; 

w 

1 %,1 



12 (9x3 

w 

2, 



12 (9x, 



jl^;^, xj 



O., 



"""g V^c/3^c,/3 ^a/S^a/S/ Qrjr.2 



l-^w) ■^t/'-^g 



, /^L*^L , .R..R\ r 1 ■9^^4.1 _ 1 ^'^3.1 _ 29 (92ft ^ _ 11 9ft., I 






3 92ft, 



.1/ X , 5%,i 



"l8 9x2 v-'^^ ' 4 ^a; 

V^a/3 ^a^ ^Cil3^cil3j { 

1 9ft,, 



}(a;w,a;J 



O. 



1 ^^^4,1 1 ^^^3,1 _ 5 92 ft^^ 

144 9x4 12 9x3 24 9x2 



12 9x„ 



-}(x„,xJCs 



_i_( f-L*/-R /■R*/-L \n^c( 

\J>ci^^cil3 ^a/3 ^a/3/ q\ 



1 %,i , 7 92 ft^ 
24 9x2 



16 9x„ 



+ 



}(^w,a;JO, 



+ 



(29) 



where elhpses represent those relatively unimportant corrections. Comparing the result 
in Eq.( l26l) . the contributions from the corresponding diagrams contain the additional sup- 
pressed factor m^ / K^ when both of virtual charged gauge bosons in Figll](a) are replaced 
with the charged Goldstone G"^. However, we should consider the corrections from those two 
loop diagrams in which one of virtual charged gauge bosons is replaced with the charged 
Goldstone G^ since it represents the longitudinal component of charged gauge boson in 
nonlinear R^ gauge. As the closed fermion loop is attached to virtual W^ gauge boson 
and charged Higgs simult aneously, the corresponding triangle diagrams belong to the fa- 
mous Barr-Zee type diagrams [131]. It is shown [iJ] that this type diagrams contribute to 
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important corrections to the effective Lagrangian. For the reason mentioned above, we also 
generahze the resuh directly to the diagrams in which a closed heavy loop is attached to 
the virtual H^ and W^ fields simultaneously. 



B. The corrections from the diagrams where a closed heavy fermion loop is at- 
tached to the virtual W^, G^ {H^) bosons 

Similarly, the renormalizable interaction among the EW charged Goldstone/Higgs 
G^ {H"^) and the heavy fermions F^^^ can be expressed in a more universal form as 



C 



S±FF 



■,R 



G-F^{g^j^u^ + g:fu:^)Fp + H-F^{n'2^uj_ + ?/:>+)F^ + h.c. , (30) 



where the concrete expressions of G'^g'^i "^^3'^ depend on the models employed in our 
calculation, the conservation of electric charge requires Q/3—Qa = 1- Generally, the couplings 
among the charged Goldstone/Higgs and quarks are written as 



C 



S±dti 



eV* ( - 
V2sJ 



m m 



u + H'd 



m 



m.. 



m m 



u] + h.c., (31) 



where the parameter Be also depends on the concrete models adopted in our analysis. In 
full theory, the couplings in Eq. fl30|) induce the corrections to the effective Lagrangian for 
6 — i- 57 through the diagrams presented in Fig{T]^b, c). 

Since there is no mixing between the charged gauge boson and charged Higgs/Goldstone 
at tree level, the corresponding corrections from the diagrams presented in FigiU^b, c) to the 
bare effective Lagrangian do not include the ultraviolet divergence, and can be formulated 
as 

^WH = Tis'^O ^'" ^"' 1 '3 ^" ^^ ' ^"^ ' ^' ' ^^" ' ^^^ 

+ (x„x,^)V2p3(x„,x^„x„a;,^,x,^)[5R(?^-X-^lX1.)^^ 

Hx^xj'/'p,{x^,x^,,x,,x,^,x,^)[^{n;'::C-n;'Xjo, 
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Hx,^xJ'/'P,{x^,x^,,x,,x,^,x,^)[^{h;X,-^X,)^s 






WG WH^ '^ ^ ^ 13a 13a ' H± ™ 



(32) 



The expressions of form factors Pi(x, y, z,u,w) {i = 1, ■ ■ ■ , 4) can be found in appendix. 

Using the asymptotic expressions of $(x, y, z) at the hmit x, y ^ z m Eq j27l we simphfy 
the expressions of Eg. (1321) in the hmit nip = irip 



nip > m„ as: 



reff ^ ^G pa Jim ^ 



r21 



+ (- + ^)fln 



X 



+ 



X 



16 48 

'C,L/-L 



m 



64 288^'^ 



m? — m^ 



f± 



IC,R/-R 



\ pa ^ap pa ctp J ^ \ pa ^aj3 pa ^ ap J -'--'- 

19 - 20g/3 ^ 2 - 4Q;3 



(in 



.2 ™2 



^2,i("^„,"^J-^2,i("^ 



m 






r/i' 



771^ 



H± 



144 48 



pl6 ^2 + 6g^. 2 

^ mm 

144 48 V F 



m,^ 



m^ 



X 






L 

Q/3 



^^fC)0. 



iO" 



lc,R/-R 



(in 



^2,l("^w'"^?)-^2,l("^L."^?) 






X 



\2Qii 6 - 2(5^ 
T44+^^'^""^" 

"^{^tC - ^I'^O^s - ^^(^^l'C - ^;f O^ii 



m? — m? 

w u± 



8^2 



1 + In m — 



2 ^2,ii"^:,"^;j-^2.ii"^;±: 



m 



771^ 



771^ 



X 



3f^(-Hf c"-, + •H^A'l)c?3 - i^nfc\ + ^'/C'llc?, 

\ pa ^ctp pet ^ ap J " \ pa ^ap pa ^apj ^ 



16 



8^2 



1 + Inm 



^2,1 (ml , m2) - ^ (m2 , m^) 



m^ — m^ 



X 



V /3a a/9 ISa^apJ « V /3a a/9 fia^aPJ 1 



4^2 



1 + In m ^ 



^2,l("^^"^t^)-^2.l("^^±,"^D 



TTl^ — 77^2 



X 



V P<^ all Pa ^Q/3/ 8 V /3a ^Q,/3 ,9a ^0/3/ 1 



AV2 



1 + \nm 



^2,l('^w'"^?)-^2,l('^^±,"^t^) 






X 



\ 13a ^al3 13a ^aP J » \ 13a ^al3 fia^afij ^■i 



(33) 



The results indicate that the corrections to the effective Lagrangian from the diagrams 
presented in Fig{T]^b, c) are suppressed in the limit m^ = nip = nip ^ m„ unless the 
couplings 'H'^'^'^ violate the decoupling theorem. 



It is wel 
NLO level 



known that the short distance QCD affects the rare B decay strongly. At the 



15| . the Wilson coefficients at the bottom quark scale are given as 



C'8M^0.7(C'8(/iJ + 0.12), 



0.88) 



(34) 



where the corresponding Wilson coefficients at EW scale are written as 

C'sl/uJ = C2(/iJ + C5(/iJ + C9(;u„) + Cn(/i„) , 

Csl/Xj = C6(/iJ + C8(/iJ + Ci2(/iJ + C^13(/iJ . 



(35) 



As an application, we investigate the relative corrections to the branching ratio of rare 
decay i? — ?■ X^7 originating from those sectors. 

III. THE CORRECTIONS TO BRANCHING RATIO OF 5 ^ X^7 

In order to eliminate the strong dependence on the b-quark mass, the branching ratio 
is usually normalized by the decay rate of the B meson semileptonic decay: 
r(i?-^X^7) r(6-s7) 2a^ 



r{B^X^eu) Tib-ceu) 3np{y)x{y) 



|C5(/i.)h 



(36) 



where p{y) = 1 
and xiy) = 1 - 



iy 



ly^ — y^ — 12y^lny is the phase-space factor with y = {m^/m^y, 



'-^^f{y) with f (ml/ml) 
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2.4. From now on we shall assume the 




(b) ^6„/7r 

FIG. 3: The relative correction to the branching ratio of the inclusive B — )• X^j decay versus the 
possible CP violation phases 0Qp. Where the solid-line represents the theoretical correction with 
Ql3 = 1, Qa = 0, the dash-line represents the theoretical correction with Q^ = —1/3, Qa = —4/3, 
the dot-line represents the theoretical correction with Q^ = 2/3, Qa = —1/3, the dash-dot- 
line represents the theoretical correction with Qjs = 4/3, Qa = 1/3, and the dash-dot-dot-line 
represents the theoretical correction with Qj^ = 5/3, Qa = 2/3, respectively. The enhancing factor 
is chosen as the trivial ;Sc = 1 in (a), or a nontrivial Be = 10 in (b) 



value BR{B — )■ X^eu) = 10.5% for the semileptonic branching ratio, a^{mj = 0.118, 
a^{m^) = 1/127. For the mass spectrum of SM, we take m^ = 174 GeV, m^ = 
4.2 GeV, m„ = 80.42 GeV and m, = 91.19 GeV. In the CKM matrix, we apply the 
Wolfenstein parameterization and set A = 0.85, A = 0.22, p = 0.22, rj = 0.35 IGj. 

Without loss of generality, we adopt the universal assumptions on those couplings and 
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mass spectrum of new physics as 

^ctp ^ap poL poL pa. pa. ' 

"^^a = "^F^ = "^«± = Kp (37) 

To continue our discussion, we assume the electric charge of heavy fermions as Q/s = 
2/3, 1, —1/3, 4/3, 5/3, which corresponds to the electric charge of another heavy fermion 
in inner loop Qa = —1/3, 0, 4/3, 1/3, 2/3 respectively. In addition, we also assume that 
those heavy fermions with fractional electric charge all take part in strong interaction. 

Many extensions of the SM include the heavy fermion fields with Qjs = 2/3, Qa = —1/3. 



In the extensions of SM with large [18| or warped [19!] extra dimensions, the KK excitations 
of up- and down-type quarks form a closed fermion loop which can be attached to the zero 
modes of charged gauge boson and Higgs. In the minimal supersymmetric extension of 



SM (MSSM) 20|], the closed fermion loop composed by chargino (Q^g = 1) and neutralino 



Qp , = 0) can be attached to the charged gauge boson and Higgs. In the 3 — 3 — 1 model 



2l| . the electric charge of exotic quarks are assigned as Q^ = 4/3, 5/3. 

In many EW extensions of the SM, the couplings among the charged Higgs and quarks 
contain an enhancing factor Be- For example. Be = tan/3 in the MSSM is a strong enhancing 



factor at large tan/3 limit. In other EW theories such as the littlest Higgs 22|, 3—3 — 1 model, 
the couplings among the charged Higgs and quarks also contain a nontrivial enhancing factor 
Bc^ I. In our numerical discussion, we assume the possible enhancing factor with a trivial 
value i3c = 1 or a nontrivial value Be = 10. 

Including NLO QCD effects, we plot the relative corrections to one loop SM theoretical 
prediction on branching ratio of the inclusive B -^ X^'-f decay versus the possible CP 
violation phase (j)^p with i3c = 1 in FIG.|3]^a). Depending on concrete choices of Qjs and the 
CP violation phase 0pp, the relative corrections to the branching ratio from those two loop 
diagrams can reach 2.5%. Comparing with the corrections from QCD, the modifications 
from those two loop EW diagrams are unimportant certainly. Nevertheless, those effects 
can be observed possibly in the experiment along with improving of the theoretical analysis 
and increasing of the experiment precision. Taking the enhancing factor Be = 10, we plot 
the relative corrections to one loop SM theoretical prediction on branching ratio of the 
inclusive B — )■ X^7 decay versus the possible CP violation phase (p^p in FIG. Et^b). Because 
the contributions from two loop Bar-Zee diagrams are enhanced drastically, the relative 
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800 1 OOO 

A,,/GeV 



FIG. 4: The relative correction to the branching ratio of the inclusive B — )• X^j decay versus the 
energy scale of new physics Aj^p. Where the solid-line represents the theoretical correction with 
QjS = 1, Qa = 0, the dash-line represents the theoretical correction with Q^ = —1/3, Qa = —4/3, 
the dot-line represents the theoretical correction with Qp = 2/3, Qa = —1/3, the dash-dot- 
line represents the theoretical correction with Q^ = 4/3, Qa = 1/3, and the dash-dot-dot-line 
represents the theoretical correction with Qp = 5/3, Qa = 2/3, respectively. 

corrections to one loop SM theoretical prediction on the branching ratio of S — ;■ X^7 can 
reach 4.5%. Although the two-loop EW corrections can not compete with that from QCD, 
we cannot neglecte the corrections with this magnitude. 

Assuming Be = 10 and ^^^p = 7r/2, we plot the relative corrections to one loop SM 
theoretical prediction on the branching ratio of 5 — )■ X^7 varying with the energy scale of 
new physics Aj^p in FIG. HI Since the intervention between the top quark and the particles in 
new physics, the relative corrections reach the maximum (~ 4.5%) around A^^p = 200 GeV. 
With increasing of A^^p, the relative corrections turn smaller and smaller. At Aj^p = 1 TeV, 
the relative corrections are about 2%. 

In the SM, the CP asymmetry of the 5 — )■ X^7 process is calculated to be rather small: 
A^p ~ 0.5% [17]. Certainly, the new CP violation phases may induce the observable effects 
on the CP asymmetry of i? — t- X^'-f. However, the numerical results indicate that the 
corrections from those two loop diagrams to the CP asymmetry of -B — ^ X^'y are rather 
small. The relative correction to one loop SM theoretical prediction on the branching ratio 



of i? — 7- X^7 is already above 8% when Be = 30, Aj^p 



200 GeV and 



'CP 



7r/2, the 



corresponding correction from those two loop diagrams to the CP asymmetry is still smaller 
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than 1% under our universal assumptions on the parameter space. 

As mentioned above, the universal assumptions on the couplings and mass spectrum of 
new physics are adopted in our numerical analysis. In concrete EW extensions of the SM, 
this choice is a very simple assumption on parameter space. However, the numerical results 
given above reflect the typical magnitude of corrections from those two loop diagrams to the 
branching ratio of 5 — )■ X^7 unless there is contingent cancelation among different sectors 
of those two loop diagrams in concrete extensions of the SM. 

IV. CONCLUSIONS 

Applying effective Lagrangian method and on-shell scheme, we analyze the EW cor- 
rections to the rare decay 6 — ?■ s + 7 from some special two loop diagrams in which a closed 
heavy fermion loop is attached to the virtual charged gauge bosons or Higgs. The analysis 
shows that the final results satisfy the decoupling theorem explicitly when the interactions 
among Higgs and heavy fermions do not contain the nondecoupling couplings. Adopting the 
universal assumptions on the relevant couplings and masses of new physics, we present the 
relative corrections from those two loop diagrams to one loop SM theoretical prediction on 
the branching ratio of 5 — )■ X^7 varying with the possible CP violation phases and energy 
scale of new physics. The numerical results indicate that the relative corrections from those 
two loop diagrams can reach 5% if there is not contingent cancelation among different sectors 
of corresponding contributions. 
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Appendix A: Form factors in the two-loop Wilson coefficients 
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The definition of "^{x, y, z) is written as: 
• A^ > 0, yy + v/i< yi: 



^!{x,y,z) = 21n 



x + y-^-X.^x_-y + z-X._^^y_^^z_ 
2x ^ ^ 2x ' XX 



-^L..{ 



X + y — z — X 



2x 



^Li,{ 



X — y + z — X\ Ti^ 

+ - 

2x ) 3 



^ , (Al) 



where Li^{x) is the spence function; 
A2 > 0, v^+ v^< ^: 



*(x,i/,z) = Eq.(|XI])(x^y) ; 



(A2) 



A^ > 0, yi+^< yi: 



*(a;,i/,z) = Eq.(inD(a;^2;) ; 



A2<0: 



\l/(x,|/,2;) = 2|C/2(2arccos 
+(7/2(2 arccos( 



-x + y + z,^ / x-y\z. 
+ Cf2(2arccos( ^ — 

^ 2^xz ' 



X \ y — z 
2Jxyy 



))} 



where C/ofx) denotes the Clausen function. 



The expressions of ipo{x,y), ipi{x,y), ip2{x,y) and ip^{x,y) are given as 



^o{x,y) 



{x + y)\Tix\n.y + {x — y)Q{x,y) , x > y ; 
2a; In^ a; , x = y ; 

{x + y)\nx\ny + (y — x)Q(y,x) , x < y . 



(A3) 



(A4) 



(A5) 



^i{x,y) 



— In X In y — |3^9(a;, y) , x > y ; 



4 — 2 In X — In x 



x = y- 



— Inxln?/ — ^3^9(?/,x) , x<y 



(A6) 



^2{x,y) 



(2x'^+&xy) In x-{6xy+2y^) In y 
{x-y)^ 

{2x^+6xy) In a:-(6xy+2y^ ) In y 
(x-y)^ 



4xy 



(x-y)'- 



■Q{x,y) , x> y 



x = y 



ixy 



{y- 



'^e{y,x) , X <y 



(A7) 
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with 



v^{x,y) 



12xy{x+y) 
(x-j/)5 



Qix,y) 



2{x'^ +6xy+y'^) 
{x-yY 



. 2(x'^+20x^y+lla:y^) \nx-2{y^ +20xy'^ + llx'^y) \ny 

(x—y)^ ' " ' 
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+ ^lnx , 



150x2 ' 5x2 

l2xy{x+y) Q/ ^\ _ 2{x^+&xy+y^) 

2{x'^+20x'^y+llxy^) lnx-2(y'^+20xy'^ + llx'^y) Iny 
{x-yY^ 



x = y; 



X <y , 



(A8) 



Q{x,y) = Inxln 21n(a; — y)ln 2Li2{—) H . 

The functions adopted in the text are written as 

x^ \yy^ X — y^ Xw^ y ^, x^^ix^u^v) — y^^iy^u^v) 
Qi,i [x, y) = , lli{x, y; u, v) = 



x-y 
Fi{x,y,z,u) = — ^2,1(2, m) 



x-y 



5^^3,1 , o-^^^a,! 



dx^ 



+ 3- 



dx^ 



{y,x) - { 



1 ^ga.i 1 '9^^3.1 

8 dz 24 dz^ 
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32 (9^2 iQ Q^3 128 9^4 
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dx^ 
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24 (9x3 



4 (9x2 
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— I (2; + -u) + 2(2;ln2; + -uln-u 
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x,y) 
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(9x3 
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(9x2 
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(9x4 
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dx^du 
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dx"^ 
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